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ABSTRACT 


A FORTRAN program has been developed to allow for the 
use of probabilistic design methods in the numerical optimi- 
zation process. The program was written as a set of sub- 
routines for COPES (Control Program For Engineering Synthesis). 
COPES maximizes or minimizes a numerically defined objective 
function subject to a set of inequality constraints using 
the optimization program CONMIN (A Fortran Program for Con- 
strained Function Minimization). The program developed here 
allows for the use of both the normal and lognormal distribu- 
tion models. Design examples are presented to demonstrate 
the program capabilities. User instructions are provided for 


inclusion in the COPES user's manual. 
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in “INER@ DUCTION 


In the traditional deterministic approach in engineering 
design the variables are treated as single valued numbers. 
These variables, whether dimensions, material properties, 
Hoads, etc., are actually statistical in nature. By the use 
of safety factors designers protect against these variations 
by usually over designing. A more logical approach is to 
take into account the statistical data known about each 
variable and design for a certain reliability. By consider- 
ing the statistical nature of each variable we should be 
able to better predict reliability and performance. 

Numerical optimization has proven to be a very powerful 
tool in engineering design. Virtually all design problems 
require minimization or maximization of some objective. For 
the design to be acceptable, it must also satisfy a certain 
set of specified requirements called constraints. If these 
constraints are specified as probabilities of failure the 
result would be a probabilistic design. 

The purpose of this research was to test the applica- 
bility of combining probabilistic design concepts with those 
of numerical optimization by: (1) developing a pilot computer 
code to calculate probabilities of failure, (2) incorporate 
this computer code with COPES/CONMIN in order to perform 
probabilistic optimization, and (3) test using numerous 


examples. 


Ae 





It 1s assumed throughout this discussion that the reader 
is familiar with the use of the COPES/CONMIN optimization 


program {Ref. 1], and [Ref. 2]. 
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If. NUMERICAL OPTIMIZATION 


The numerical optimization problem considered here is 
stated as follows: Find the set of n design variables 


contained in the vector X which will 


Minimize F(X) (As) 


Subject to: 


Sia 2S) ee J eae a (2) 
xX < x. < xl i = \ ten (3) 
i in- ii-r- i 


The components, Xr Oreo meweerCrued sto aS design vagiaples 
Pelchedce changed to improve the design. The function F(X) 
1s called the objective. Inequality constraints, Bia WES are 
the response limits imposed on the design. There are n de- 
Sign variables and m inequality constraints. The lower and 
upper bounds, 5) and Se are limits imposed on the design 
variables to insure a practical result. 

For a deterministic design COPES [Ref. 1], determines 
the constraints and objective as follows. Consider the 


design of a single bar undergoing uniaxial tension. If we 


wish to minimize the weight, the objective, 


F(X) = pAL (4) 


ie 





where 


0 = the specific weight; 
A = the area of the bar; 
L = the length of the bar. 


The design variable is the area, A. If the stress in the 
bar is S and the stress limit imposed is Sy we desire that 


S be less than or equal to Sy so the constraint becomes 


Tova = 20 (5) 


In this investigation a computer subroutine was developed 
to provide constraints based on probabilities of failure. 

If the allowed probability of failure is PF and the 
actual probability of failure was calculated to be Pf, then 


the constraint would be 


P£/PF - 1 < 0 (6) 
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PPh r RObAR Lio te = Dn Sua 


A. RANDOM VARIABLES 
If one has a large population of n elements for which 


some parameter X,,isl,n is determined, the mean value of 


the population is 


1 2 
ea ae = 2} (x; ) (7) 
i=l 
The variance is 
n 
_ tL 2 
is) x) (8) 
a 
The standard deviation is 
Oo = /V_ (9) 


The standard deviation 1S a measure of the amount of varia- 
bility of the population data. An additional measure of 
the variability is the dimensionless coefficient of variation. 


In engineering design the coefficient of variation is fre- 


quently used. 


i Be 





ae SeINCTIONS OF SEVERAL VARIABLES 


If Y is a function of n variables, the standard deviation 


of the function Y can be approximated as follows: 
If Y = £(x,; 1 = l1,n) 


Assuming that all xX. are independent random variables 


mien a@cording to [Ref. 3, p. 59] 


Ben 2? we 
Oy ee Ll cra, O d (11) 
ae a 1 


lo-19 


i 


on STATISTICAL MODELS 
Ll. Normal Distribution 
The normal distribution model is one of the most 


widely used. The normal density function is given by 





2 
ae tS) 
ee ee XP > (e2) 
/270 ZG 
xx x 
where: 
£y (x) =eeenemoropoDTiliityecensity Lunction of 
the random variable X; 
Oe. = standard deviation of the random 


variable X; 


u = mean of the random variable X. 


16 





The magnitude of the standard deviation determines the dis- 
persion of the distribution. A large standard deviation will 
result in a wide bell curve and a small standard deviation 
will result in a narrow one. The curve is symmetric about 
the mean. Approximately 68.26% of all samples will lie 
within plus or minus one standard deviation of the mean value. 
The standardized form of the normal density function [Ref. 4: 


p. 194], can be expressed as 





1 wee 2 
ON ee Same exp [——] (3) 
¥2n 
where 

x - WU 
eo (14) 

x 

o, (y) = the standardized form of the single 


variate normal density function. 


The single variate cumulative normal distribution function 


1S given by 


Xx 
= 5 
Ba) J £,.(x) dx (BS) 
where: 
F(x) = the probability that the random variable 
as X is equal to or less than the specific 


value xX. 


Ig 





Or, uSing the expression, 





Equation (15) can be rewritten as, [Ref. 4: p. 195] 





1 (x-n)/o,, 2 
a | exp ([—4—] dy (16) 
¥27 -% 
or 
Ol 
_ Ded 
F(x) = ol 5, ] (17) 


2. Lognormal Distribution 


Tf ln(X) has a normal distribution then X is said to 
have a lognormal distribution. The lognormal density func- 


mEeneeimet. 4: p. 196], is given by 





ee 
f(x) = —S— 2 exp(-nt ), xo (18) 
J27c.. * 2 
x x 


where: 


mean of In(X); 


SS) 
i 


> 


Seaneard deviation o£ In(X). 


18 





The cumulative lognormal distribution PUMCEVOn, [ReE. 4: 


pe 196], is given by 


a 


oe ee. 
no 1 f ox -y? 
F x = —- exp [—-—] dy; xe 0 1 
x V2 -0 ss 
Or 

_ in (x) -X 

Fy (x) = oe a: mae 20) 
“ 


The mean, eee and the variance Vt of X are 


ce 
Bae allen exp[X + =!) ZA) 
V = eee iexo (6. ould (22)) 
x X x 
solving for X and o 
o. i 2 
a in(yy) Zz in (Vy + uy) (Zs) 
6, = ,/-1n(u2) + In{v, + u71 (24) 
x x K se 


D. PROBABILITY OF FAILURE 
Consider the constant area bar in Figure 3.1 under load 


P. If the stress is the bar is S and the stress limit in 


ig 





Figure 3.1. BAR UNDER UNIAXIAL LOAD 


tension for the material is Sy, then the probability of 
failure of the bar is the probability that the stress S$ will 
be greater than the limit Sy. The method for determining 
the probability of failure in the bar for both the normal 
and the lognormal formats 1S presented below. 
De Netanal Distribution 

If both S and Sy are random variables which follow a 
normal distribution then the probability of failure is the 
probability the S$ is greater than Sy. The failure function, 


[Ref. 5: p. 8-24] 


nS Ss) (25) 


20 





is also normally distributed. The probability of failure 


is the probability that the failure function is less than 


Se ecual to 0.0 or 


Ee (oa 
In the standard form 
: ees 
Baa Calas ] 
Zz 


Miemevemy the algebra of functions, [Ref. 3] 


and 


Mtewoaarery index, [Ref. 5: p. 8-24], is defined as 


BMicwewoedbality Of failure is then 


a == ?,[-8] 


2 


(26) 


(27) 


(28) 


(22) 


(SO) 


cei), 





Seeme abea under the unit@mormal density function from 6 


ZO CO. 


co il _~4 
2) xD | oe (22) 
B V2 


Zeer Lognermal Distribution 
If the same two variables S and Sy follow a logri- 
thetic normal distribution then the probability of failure 


1s expressed by, [Ref. 5: p. 8-27] 


Pie ee <5) = 8Ptintsy/s) <0) (e23) 
Letting 
Z = In(Sy/S) (34) 
Gaon 
Pf = P(Z < 0) (33) 
Or 
P£ = 1-8] (36) 


Mpeometecty index, (Ref. 5: p. 8-29] is 


Z2 





Be iIn(Sy/S) (37) 


2 2 
ah Ce (1+ Co) ] 


where the Sy and § are median values and 


i. u 
2 
ee CSy 
2 u 
1 +c? 
S 
where 
“sy 
= (40) 
SY Usy 
O 
ee < = (41) 
S 


5 will be normally distributed and the probability of 


taidware is 


oo o, | (42) 
i 2 


Ih eG lias6 
Oi I Sy) | ail 
The probability that S is greater than Sy is then determined 


in the same manner as the normally distributed case, by 


integrating the unit normal density function from 8 to ~. 
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IV. PROBABILISTIC OPTIMIZATION 


HiemOsvecciva Of probabilistic optimization is to mini- 
mize or maximize the mean value of a given function subject 
to a set of constraints based on allowed probabilities of 


failure. 


mee CONSTRAINTS 

The computer subroutine developed in this investigation, 
COPE19, calculates the constraints as described in Chapter 
II. Given a constrained variable in terms of a desired 
probability of failure, one needs only to calculate the 
actual probability of failure to determine the constraint 
value. 

Consider the cantilevered beam in Figure 4.1. A proba- 
bility of failure, PF, is desired for stress in bending. 
The mean values and coefficients of variation for the yield 
strength, Sy, the length L, the load P, and the dimensions B 


and H are given. The stress in bending 


Soa. (43) 


The standard deviation of the bending stress 


a Sle inssh 2 0SE. 2 
oSB 44 
ce (—SOp) + (aE) + (ap oy) + (FRB) (44) 
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Figure 4.1. CANTILEVERED BEAM 


z 


> 8 


as 


Zs 





The safety index 


Piescalembated probability of failure 


Pf = P, 1-8] (46) 


Bier cONnsSEraint will be 


In the computer program developed, COPE19, the required 
partial derivatives of the constrained variables with 
respect to the variables in the expression are calculated 


by finite difference methods. 


B. OPTIMIZATION 

Consider the three bar truss in Figure 4.2. The geometry 
1s specified and there are two independent load conditions. 
The design task is to determine the areas required for each 
bar that will yield a minimum mean value of the structure 
weight while keeping the stress in any member less than the 


yield stress. 
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THREE BAR TRUSS 


he 
\ WY WY 


Al AZ NS 


po P J 


Figure 4.2. THREE BAR TRUSS 


out 





Minimize: W = pH(V2Al + A2 + /2A3) 
Subject to: Sy compression S eles < Sy tension 


where S oneiae is the stress in member i under load condition We 


Pewee veaterial ~pee=.l ib./cou. in. 

Geometry H = 10.0 in. 

Loads Pl = P2 = 20,000 lb. 

Stress limits: Sy compression = -15,000 psi 
SVacenston “= 20,000 psi 

Maintain symmetry Al = A3 


Beginning with the design Al = A2 = A3 = 1.0 sq. in. and 
Optimizing using COPES/CONMIN (deterministic) the following 


results are obtained: 


We = 2,632 1b. 
Al = A3 = .7796 sq. in. 
Ae 2427/5 SC. ain 


The critical constraints are SIG11 and SIG32, which are 
equal. The design space is shown in Figure 4.3. As no 


safety factor was used, the above deterministic design 


Zs 





IN. 2} 
gi 
& 
Be 
“gt -— 
oe 
@ 


Ae (SQ. 





3 J 
mi €S@. IN, 3 


COUSTANT OBJ ec elivE. (LES) 


meee a =—— «8 CONSTANT STRESS SIGL1 (PSI) 


Figure 4.3. THREE BAR TRUSS DESIGN SPACE (DETERMINISTIC) 
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corresponds to a probability of failure of 0.5 or a 503 
chance that SIG11 and SIG32 will be greater than the stress 
limit Sy. Utilizing the subroutine COPE19 to provide proba- 
bilistic based constraints and assuming that the coefficients 
of variation of all the variables are 0.1, the following 
results are obtained when designing for a maximum probability 


of failure of 0.5 in each member; assuming a normal distribution, 


Wo] 27.6 Sen 1b. 
P= AS 7 See sq. in. 
AZ = .4094 sq. in. 


The design space is shown in Figure 4.4. The lines of con- 
stant probability of failure for SIGl1 are shown. As one 
would expect the line of constant 0.5 probability of failure 
for SIGl1 in Figure 4.4 is the same as the constant 20,000 
psi line for SIGll in Figure 4.3. The reverse is also true. 
If the three bar truss is optimized for a maximum probability 
Stfierarlure of 0.4 the resulting safety factor for SIGl1 is 
ive Oonel,501. psi and the minimum weight is 2.726 lb. Using 
19,361 psi as the limit of stress for a deterministic optimi- 
zation yields approxiamtely the same minimum weight. In both 
cases the feasible design space is the same with the line 
Seecedseant 0.4 probability of failure and constant 19,361 


PoeeelvinG On top of each other. 
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Figure 4.4. THREE BAR TRUSS DESIGN SPACE (PROBABILISTIC) 
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Table I shows the results of minimizing the weight of 
Peiewtnree bar truss for various probabilities of failure. 
Table II provides the deterministic optimum designs for the 
equivalent factor of safety. 


The difference in designs is 


within the numerical accuracy of the optimization. 


TABLE I 


Tiree BAR TRUSS PROBABILISTIC DESIGNS 


PROBABILITY WEIGHT Al = A3 A2 SiGiak SAFETY 
OF FAILURE 1S SOc INe SC ALN. leroy FACTOR 
aes 20 0 re 2) - 40938 20000 Lae0 
4 2477260 FoLoUsS ~-41785 T9861 8) 
ee Dh OM OVE . 84478 -43438 Foie ie) 7 
we 2 S30 .86887 -48544 17940 eee 
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It should be noted that the correspondence of the proba- 


bility of failure and the safety factor shown here is unique 
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eoeeiie Problem, in which all C = 0.1. In Gevuewmaie Nome inieat 


Seermrzation, this one to one relationship cannot be assured. 


C. EFFECTS OF COEFFICIENTS OF VARIATION ON THE DESIGN SPACE 

As the coefficients of variation of the variables used in 
the design subroutine become smaller, the bands in the design 
Space for the constrained values between Pf = 1.0 and Pf = 0.0 
becomes narrower. Figure 4.5 shows the design space for the 
PreeemeameeruSS in the vicinity of the optimum for the case 
where all variables have coefficients of variation, C = 0.1. 
Figure 4.6 shows the same space for the case where all coeffi- 
€wenes Of Variation, C = 0.001. In both figures the lines of 
SComsecant Goyective and probabilities of failure for the 
critical constrained variable, SIGl11, are shown. The line 
of constant Pf = 0.5 is the same for both but the band between 
PE = 0.9 and Pf = 0.1 is much narrower for the case where 
eo = 0.000: 

This phenomenon can cause numerical difficulties in the 
optimization process. Premature termination of optimization 
Can occur as the optimizer may be unable to move down the 
resulting 'narrow valley'. The absolute and relative termina- 
tion criteria of CONMIN may have been satisfied but the true 
optimum might not have been reached. 

imenie initial design 1s a considerable distance from 
the narrow band of changing probabilities of failure, CONMIN 


may obtain zero gradients for the constraints which would 
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also cause termination of the optimization process. A 
routine to attempt to prevent this was included in COPE]9. 
The cumulative distribution curve for the normal format is 
shown in Figure 4.7. The numerical integration routine used 
to calculate the area under the unit normal distribution 
Curve yields a probability of failure, Pf of 1.0000000 for 
eoeormm—e «2277921 and a Pisot 0.0000000 for a 8 of 5.5000086. 
It was modified to "widen" the zone between Pf = 1.0 and 
Pf = 0. A polynomial (a line equation) was placed at each 
ereason tne cumulative distribution routine to yield a Pf 
of 1.0 at 8 equals -25.0, decreasing to a Pf of 0.9999866 
at -4.2 and a Pf of 0.00001335144 at 8 equals 4.2 decreasing 
Bema Ore0.0 at 25.0. This is shown in Figure 4.8. This 
provides the same optimization results within the accuracy 
of the program, but avoids the problem of zero gradients of 
the constrained variables. 

An additional approximation was used to determine the 
probability of failure when the coefficients of variation 
concerned are zero or nearly so. The probability that xX is 


greater than Y is, in the standard format 
u 
oie a Oe] (48) 


where 


(49) 
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Figure 4.8. MODIFIED CUMULATIVE NORMAL DISTRIBUTION CURVE 
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O = Oo +06 Call) 


Mathematically, if a — 0 7) Een One (On thnee scituataons (can 


eccur. 
mo. —- Y, 6B = 0, then Pf = 0.5; 
mmeexX et, 6 > -—>, then Pf = 1.0: 
ee = Yt, PB > ©, then Pf = 0:0. 


@aas 2s shown graphically in Figure 4.9. In order to provide 
CONMIN with a smooth function with which to work, a polynomial 
was used to connect the Pf = 1.0 and Pf = 0.0 in Figure 4.9. 
This 1s shown in Figure 4.10. The interpolating polynomial 

is used when 8, the modified safety index, is within one 


unit of 0. Where 8 is defined to be 


CT ae Sioa? Soe (51) 


A similar approximation was developed for the lognormal 


model. 
Experience has shown that these modifications dramatically 
improve the numerical stability of the optimization process 


without significantly affecting the mathematical result. 
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V. COMPUTER PROGRAM 


A. COPES INPUT DATA 

Standard deterministic optimization requires that COPES 
data blocks A-0 and V be used. For probabilistic optimiza- 
tion the same data is entered in these blocks with the 
exception of NCALC and IPROB in data block B. Here an NCALC 
ee eveula be (entered to indicate probabilistic optimization 
was to be performed and an IPROB = 1 or 2 to indicate which 
probability model is to be used. Additional data blocks are 
required for probabilistic optimization. Appendix A con- 
tains detailed input instructions for these data blocks as 


well as an example problem including sample input data. 


B. OPTIMIZATION DATA FLOW 

The interaction between COPES, CONMIN, and the analysis 
subroutine will be the same for probabilistic optimization 
as it is in standard deterministic optimization with one 
exception. In the optimization process, each time COPES 
requires the determination of a constraint vector it will 
call COPE 19. COPE19 will call the analysis subroutine as 
many times as necessary to determine the required partial 
derivatives by finite difference steps. Once these are 
obtained, the constraint vector is determined and provided 
to COPES. A flow diagram for probabilistic optimization 


is provided in Figure 5.1. 
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Figure 5.1. PROBABILISTIC OPTIMIZATION FLOW DIAGRAM 








S- ©OrlEt9 SUBROUTINE 


The main computer subroutine developed in this investi- 


gation is COPE19. Three additional subroutines, COPE20 


through COPE22, were also developed to perform specific 


tasks required by COPE19. The operation and flow of these 


subroutines are discussed below. 


TP eput 
The following is provided to COPEL9: 


a. 


DD. 


The COPES control parameter ICALC. 

Which probability model to use, either normal 
or lognormal. 

The number of probability variables which make 
up Ehe constraint equations (IVAR). 

The global locations in the ANALIZ subroutine 
common block of each IVAR. 

Coefficient of variation of each IVAR, assumed 
Gousrtam: LunTOugGnouc the Optimization. 

The number of constrained variables. 

The global location in the ANALIZ subroutine 
common block of each of the constrained 
variables. 

The upper and lower limits imposed on each of 
the constrained variables. 

The allowed probability of failure at the upper 
and lower limits of each constrained variable. 
The coefficient of variation of the constrained 


variable limits. 
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K. The variable MGRAD which determines how often 
the constrained variable gradients are calculated. 
2. Determination of Constraint Vector 
When COPES requires a constraint vector for a given 
design it will call COPE19 with ICALC = 2. COPEL9 will then 
calculate the partial derivatives of each constrained varia- 
ble with respect to the variables which make up the con- 
straint equations. The subroutine COPE20 controls the finite 
difference calculations. In the example in Section ny 22, 


the cantilevered beam, it would calculate: 








ISB 6P : 
oS (52) 
BH 
ISB ay 
BH 
ISB L2PL 
i eae 9a: Jo) 
BH 
3SB 6PL 
a ee (35) 
B“H 


COPELIY will then calculate the standard deviation (Equation 
44) of the constrained variable, in this case SB, as dis- 
cussed in Section III.B and Section IV.A. Now the safety 
index is calculated and subroutine COPE21 called to deter- 
mine the probability of failure. When required, COPE21 will 
calimesiotouLtine COPE22 to perform, numerical integration of 


the normal density function. The constraint value to be 
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stored in the G vector is then determined (Equation 6). 
The above would be performed for each constrained variable 


and the resulting G vector provided to COPES. 


D. REDUCTION OF FUNCTION EVALUATIONS 

In order to determine the gradients of the constrained 
variables finite difference steps are taken resulting in a 
Significant increase in the number of analyses performed in 
the optimization process. In many problems these gradients 
remain essentially the same throughout the optimization. 
Therefore, it is sometimes possible to reduce the frequency 
with which the gradients are calculated and achieve approxi- 
mately the same result. The input variable code LGRAD is 
included as an option in probabilistic optimization to allow 


this choice. LGRAD has the following meanings: 


0: Calculate gradients each time COPE19 is called. 


1: Calculate gradients at the beginning of each 
CONMIN iteration. 


Peco leuilate gradients only at the beginning OF 
Seem aec von - 


Table III demonstrates the result of various LGRAD selections 
when designing the three bar truss for an allowed probability 
of failure of 0.50, where all coefficients of variations, 
a=. Ose. 

Probabilistic optimization was performed using the three 


LGRAD values on numerous additional test cases. In almost 
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TABLE Ii 


REDUCTION OF FUNCTION EVALUATIONS: THREE BAR TRUSS DESIGN 


LGRAD WEIGHT Pas = a3 PZ NUMBER OF 
LBS. SO. mi. S@. HN. FUNCT EVALS 
Deterministic 2.6326 .77962 82752 35 
Prebabdilistic 0 2.6390 78827 40938 343 
Probabilistic 1 2.63189 . 77386 25011 244 
Pxrebabilistic 2 2.035 .78972 s4042 7 bs. 


every case the optimization results were essentially the same 
for LGRAD = 0 anc 1. The composite driveshaft design in 
Giapter V provides an excellent example showing the fumction 
evaluation reductions when performing a relatively complex 
Gesiagn. Table IV demonstrates the results of reducing the 
frequency c= constrained variable gradient calculations 

wnen Gesicning the composite driveshaft for minimum weicht, 
Meaac cle nofmal distribution model. For this particular 
probiem, calculation of comstrained variable gradients only 
at the beginning of optimization was insuificient to cbtain 


results near the ootimun. 


TABLE IV 


REDUCTION OF FUNCTION EVALUATIONS: DRIVESHAFT DESIGN 


LGRAD WEIGHT NUMBER OF FUNCTION 
BES . EVALUATIONS 
0 97,0203 TOME 
il S012 05 1344 
2 Zac LG 89 
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VI. NUMERICAL EXAMPLES 


The following examples are presented to demonstrate the 
capabilities of the computer subroutines developed in this 
investigation. In both cases an existing ANALIZ subroutine 
was used after minor modifications. For these two shaft 
designs the subroutine from [Ref. 6], was used. The ANALIZ 
Subroutine was altered to include the required probability 
variables in the global common block and to remove the 
factor of safety calculations which were originally used as 
constraints. 

In each example a similar deterministic design was per- 
formed to give the reader a better feel for the results. 
The deterministic designs are not intended to duplicate the 
Eaovdeitlistic design results. 

The mean values used in the designs are from [Ref. 6]. 
The coefficients of variation are reasonable assumed values 


used for demonstration purposes. 


am. STEEL DRIVESHAFT 

Design, for minimum weight, a steel driveshaft to trans- 
mit 150.0 horsepower at 300.0 RPM. The loading is presented 
in Figure 6.1. The shaft is designed against failure se 
strain, in torsional buckling, vibration frequency anda 
maximum deflection of 0.05 inches. 

The design variables are the shaft thickness and the 


shaft inside diameter. The thickness is to be maintained 
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between 0.01 and 2.0 inches and the inside diameter between 
Q0.5 inches and 5.0 inches. The initial design is an inside 
diameter of 2.0 inches and a thickness of 1.0 inen. Pie 


remaining data for the analysis is presented in Table v. 


TABLE V 


INPUT DATA FOR STEEL DRIVESHAFT 


MEAN COEP Se veh r 
VALUE OF VARIATION 
HORSEPOWER SOT US 
SPEED 300. RPM S16) 
AXIAL LOAD 50. EB. eS 
LENGTH 120 in ks 
YOUNGS MODULUS 30,000,000 PSI .038 
SHEAR MODULUS io oO O. a PST .042 
POISSONS RATIO — .026 
SoBeirrTC WEIGHT Seo eS eure iN SE 
STRAIN LIMIT =o OEE Serb) Wee aot 
MAX DEFLECTION 205 IN. Ae 
THICKNESS 2025 
INSIDE DIAMETER ag OG Ue: 
FAILURE DUE TO: 
TORSIONAL BUCKLING aad 
FREQUENCY .06 


Allowing for maximum probabilities of failure of 0.01 in 
strain, torsional buckling, deflection and vibration fre- 
quency; a minimum weight of 54.15 lbs is achieved for the 
normal model and 52.88 lbs for the lognormal model. The 
initial design summary is presented in Table VI and the 
final designs summaries and COPES optimization results in 
maotes Vit through XII. 

Designing for a factor of safety against failure of 2.0 


feeeserain, torsional buckling, vibration frequency, and 
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TABLE VI 


STEEL DRIVESHAFT INITIAL DESIGN SUMMARY 


STEEL DRIVESHAFT OUTPUT 
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OPTIMIZATION RESULTS 
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TABLE IX 


STEEL DRIVESHAFT DESIGN SUMMARY: NORMAL 


Sees DRIVESHAFT OUTPUT 


NUMBER OF PLYS 2 1 

NUMBER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

ECCENTRICITY = 3.0 

DIMENSIONS 

PLY THICKNESS PERCENT DIAMETER THETA 

INSIDE 0.50000E+01 
1 6.99883Z2-01 100.00  0.51998E+01 0.9 

STIFFNESS 

AR = 0.52757E+08 

PI = 0.171588+09 

GI = 0.12011F+09 
LOADS: 
a T . ? fe RPM 
1 0.31573E+05 0.0 0.50000B+02 9.15000E+03 0.300002+03 
ee SPT S.F BPLT SF 
ELY ew @ deed e es ia e e 
1 0.9477HF-06 100.00 0. 100.90 0.68213E-03 1.91 

CRITICAL SPEED = 0.250687+04 

SRI TTGR DBETECTION = 017309¢8202 

WEIGHT = 9254154R+02 

VOLUME = 0.192 04£+03 
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LOGNORMAL 


TABLE X 
FUNCTION VALUE 0.52882E+02 
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LOGNORMAL 


STEEL DRIVESHAFT COPES PROBABILISIIC OUTPUT 
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TABLE XII 


STEEL DRIVESHAFT DESIGN SUMMARY: LOGNORMAL 


SxobL DRIVESHAFT OUTPUT 


NUMBER OF PLYS = 1 

NUMEER OF MATERIAL TYPES = 1 

NUMBER OF LOAD CONDITIONS = 1 

BEeCey CT RECITY = 0.0 

DIMENSIONS 

ey Ete KNESS Peace Eur DIAMETER Hig S124 

INSIDE 0.50 000E+01 
iveorocis=-01. 100.00 e252 Sao Oa.0 

Sl eR NESS 

AE = 0.515185+08 

EI = 0.16/740E+09 

GJ = 0.117188+09 
BOA DS : 
isc « iB M yi H.P. RPM 

mecca spe 05> OS) USS UCOETO2 Oe Ta U00E+0S Os 2 00 UUEeus 

LOAD CONDITION 1 
Pey EPL Sf - oe 1k Se F. op P ak Sent. 
1 O.97054E-06 100.900 0.0 100. 00 Umea soso 5-05 Lieto 

Cort. Cal SPEED = Ua ous o nt OY 

MAXIMUM DEFLECTION = 0.72144E-Q2 

WEIGHT = Omer seers 0 2 

vORUME = ey et US 
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deflection an optimum of 73.78 lbs is obtained. The COPES 


Optimization results and final design summary are contained 


in Tables XIII and XIV. 


B. COMPOSITE DRIVESHAFT 

Design, for minimum weight a 120 in. long, four ply 
graphite epoxy driveshaft to transmit 150.0 horsepower at 
300 RPM. Figure 6.1 shows the shaft loading condition. 

The shaft is designed against failure in transverse, logitu- 
eeiaal, ame shear strain; torsional buckling, vibration fre- 
quency and a maximum deflection of 0.05 inches. Table XV 
Comtains the initial design. 

The design variables are the thickness of each ply, the 
Orientation of each ply and the inside diameter. Addition- 
ally, the thickness of ply 2 and ply 3 must remain equal to 
each other and the orientation of ply 3 equals the negative 
Gf ply 2. The following additional constraints are placed 
on the design variables: The inside diameter must remain 
between 0.5 inches and 5.0 inches and the thickness of any 
given ply between 0.01 and 0.5 inches. Ply 2 is allowed to 
vary between 0.0 and 90.0 degrees. The orientation of plys 
1 and 4 remain constant. The remaining input data is pre- 
sented in Table XVI. 

Allowing for maximum probabilities of failure of 0.01 in 
Semscain, torsiomal buckling, deflection, and vibration fre- 
quency; an optimum weight of 9.02 lbs is achieved using the 


normal distribution model and 10.28 lbs for the lognormal 
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TABLE XIV 


STEEL DRIVESHAFT DESIGN SUMMARY: DETERMINISTIC 


STEEL DRIVESHAFT OUTPUT 
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model. The ANALIZ subroutine design summary for the initial 
design is presented in Table XVII. Tables XVIII through 
Merit contain the design Summaries and final COPES optimiza- 
tion results. 

Designing for factors of safety against failure of 2.0 
mer Straln, torsional buckling, sacle hom, and vibration 
frequency an optimum weight of 10.23 lbs is obtained. Tables 


XXIV and XXV present the deterministic design results. 
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VII. CONCLUSIONS 


Numerical optimization using probabilistic design tech- 
niques provides an effective method for designing structures 
based on allowed component reliability. 

Existing ANALIZ subroutines can be easily modified to 
perform probabilistic optimization and as demonstrated, 
complex designs can be accomplished. 

The major drawback to the program developed is the large 
increase in function evaluations as compared to the standard 
deterministic optimization. This disadvantage can be 
partially offset by reducing the frequency in which the 


constrained variable gradients are calculated. 
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VIII. “RECOMMENDATIONS 


A method of scaling the variables should be investigated 
Pemceauce ili-conditioning in the optimization process. 

This is a possible area to pursue in order to alleviate 
problems caused by very small coefficients of variation. 

Further efforts should be undertaken to accomplish the 
reduction of function evaluations while maintaining suitable 
optimization results. The altering of the frequency in which 
constrained variable gradients are calculated appears to be 
the most promising approach. 

In order to reduce function evaluations it is recommended 
that COPES be modified to allow the user to supply mathemati- 
cal expressions for the partial derivatives of the constrained 
variables via the ANALIZ subroutine. These could be used 
instead of the finite different gradients calculated by 
COPE19/COPE20. Modifications to allow the user to supply 
precalculated gradients would also be helpful to this end. 
The sophisticated user can accomplish this now by using CONMIN 
Peace in CON JuUNnCL1ON with COPEL9, COPE2Z1, and COPE22. 

It is further recommended that additional work be under- 
taken to modify COPES to allow for the use of the following 
Options with probabilistic design: A) two variable function 
space; B) sensitivity analysis; C) optimum sensitivity. 


An additional worthwhile modification would be to provide 


ed 





for the calculation in COPE19 and subsequent output of the 
standard deviation of the objective function mean value. 
This investigation dealt only with component reliability 
in the design process. Work should be undertaken to ascer- 
tain the potential for the inclusion of system reliability 


considerations in the optimization process. 
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APPENDIX A 


COPES MANUAL ADDENDUM 


This appendix is intended as an addendum to the COPES 


Manual, [Ref. lL]. 


I. INTRODUCTION 


The purpose of this document is to provide user instruc- 
tions for performing numerical optimization using probabilis- 
tic design techniques. Four subroutines have been included 
in COPES (Control Program for Engineering Synthesis) to 
provide this additional design tool. 

This discussion describes the capabilities of probabilis- 
tic optimization using the COPES/CONMIN program. A simple 


design example is first presented to demonstrate the program 


Capabilities. Guidelines are given for writing analysis 
codes. The data organization is outlined and sample data is 
presented. 


This publication was written to serve as an addendum to 
the COPES Manual; "COPES--A Fortran Program for Engineering 
Synthesis," L. E. Madsen and G. N. Vanderplaats, NPS69-81- 
003, Naval Postgraduate School, Monterey, California, March 


OZ. 
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II. DESIGN EXAMPLE 


It 1S required to design the cantilevered beam shown in 
Figure 1. The objective is to find the minimum mean volume 
of material which will support the concentrated load and 


maintain an allowed probability of failure of 103. 


y y 
_— ax 
L=2001IN— Dak 


Figure 1. CANTILEVERED BEAM 
Maat 1S, 
Minimize volume = B* H *# L (LE) 
The mean bending stress in the beam must not exceed its 
limit which is 20,000 psi which has a coefficient of variation 


@ueO..0 7 


pstRES = “CG = SF < 20,000 (2) 
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The mean shear stress must not exceed its limit of 10,000 


psi which has a coefficient of variation of 0.05; 
SHRSTR = 32 = ZF < 10,000 (3) 


and the deflection under the load must not exceed one iaeh: 


G@eerficient of variation equals 0.02; 





Additionally, geometric limits are imposed on the mean 


dimensions so that; 


O25 < mB. < 5.20 (5) 
Os" H  < 20.0 (6) 
= L050 (7) 


The manufacturing procedure is such that all dimensions have 
a coefficient of variation of 0.01. The mean value of the 
eeac, loaqe2s 057000 Ib. with a coefficient of variation of 
0.07. The Young's modulus mean value is 30.E+06 psi with a 
OPlcrcecrterelent Of variation. 

The ANALIZ subroutine on page ll of the COPES manual 


is used for the analysis. Only one modification is needed. 
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The variable P, for the load, is added to the global common 
block. 

The COPES data used for the standard optimization, pages 
60 and 61 of the COPES manual is modified to perform proba- 
balastic optimiZation. NCALC in data block B is changed to 
7 to indicate probabilistic optimization is to be performed 
and IPROB equals l is entered in column eight so the normal 
distribution model will be used. Three additional data 


blocks are required, Ul, U2 and U3. 


DATA BLOCK UL: 

There are five probability variables which make up the 
Senrstraint equations 2, 3, 4 and 7. They are B, H, E, AL, 
and P; so IVAR equals 5. It is desired that gradients of 
the constrained variables only be calculated at the beginning 
of each CONMIN iteration, therefore LGRAD equals 1 is entered 


mn column two. 


$ DATA BLOCK Ul 


4 


DATA BLOCK U2: 
The global location of the probability variables and 
their corresponding coefficients of variation are entered 


mites OlLock. 


S$ DATA BLOCK U2 
iO 1 width, 3B 


2 weal height, H 
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iO. 05 load, P 
87.06 Young's mod,” & 


oy. OL length, AL 


DATA BLOCK U3: 

This block contains the allowed probability of failure 
and coefficient of variation at both the upper and lower 
bound of each constrained variable. All four constrained 


variables have upper bounds, but are unbounded on the lower 


end. 

S$ DATA BLOCK U3 

we O7, <1, .07 Boleks 
ies, .1,.05 SHRSTR 
imeew2,.L;.02 DELTA 

ieee, .1,.d01 H/B 


The allowed probability of failure at the upper bound and 
the coefficient of variation are listed in columns three 
and four. Because there are no lower bounds the resulting 
probability of failure will be zero. These first two colums 
will be ignored during probability of failure calculations 
om, GCOPES. 

Coupling the ANALIZ subroutine to COPES yields the 


following optimum design: 


Ge, 





CANTILEVERED BEAM 


AL = 200.0 

P = 0.10000E+05 
E = 0.30000E+08 
B = ie. 92 

H = 18.75 

VOL = 7261.54 
Boles — = iol 7oe2n +0 > 
SHRSTR = 0O.41727E+03 
DELTA = 0Q0.84382E+00 
H/B = 2 eet 


The COPES output for this example is provided in Figure 2. 


IfI. PROGRAMMING GUIDELINES 


The programming guidelines presented in Chapter II of 
the COPES manual all apply in performing probabilistic 
optimization. One additional rule must be followed when 
writing the analysis subroutine: All variables in the 
constraint equations which have probabilistic distributions 
must be included in the global common block of the analysis 


subroutine. 


IV. DATA BLOCKS 


iDOnen'niseauertons fOr the probabitistic optimization 


data blocks are presented in the following five pages. 
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V. NOTES TO THE USER 


Peele probabilistic optimization feature is not 
compatible with the following COPES options: 
A. sensitivity analysis 
B. two variable function space 
C. optimum sensitivity 
D. approximate optimization 

2. It is not recommended that one attempt to design for 
allowed probabilities of failure of less than 0.00001 or 
Meeaeet enan .99999 due to modifications made to "widen" the 
cumulative normal distribution curve. 

3. The optimization process is very sensitive to the 
coefficients of variation of the probability variables which 
make up the constraint equations. Coefficients of variation 
of less than 0.0001 may result in erroneous results if the 
coefficients of variation of some of the constraint limits 
are also less than 0.0001. 

4. In some cases the initial design will be such that 
the optimization process terminates unsatisfactorily due to 
he calculation of zero constraint gradients. Termination 
well occur either due to an inability to achieve a feasible 
design or CONMIN will reach the maximum number of iterations 
no matter how large ITMAX is made. This probably will be 
due to the fact that the probabilities for failure calcu- 
lated are either 1.0 or 0.0 and the finite difference steps 


are insufficient to reach the region along the cumulative 
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normal distribution curve between Pf = 1.0 and Pf = 0.0. 
This could be caused by an initial design which is a great 
distance from the optimum or by very small GGectrrcients OF 
variation. If this happens, perform a deterministic optimi- 
zation to achieve a design in the vicinity in the design 
Space desired and use this result for the initial design in 
Beevasilistic optimization. 

5. In order to determine probabilities of failure, COPES 
calculates the gradients of the constrained variables with 
respect to the probability variables which make up the 
constraint equations. This is done by finite difference 
steps resulting in numerous calls to the ANALIZ subroutine. 
In many problems these gradients remain essentially the same 
so the user is given the ability to determine the frequency 
in which these gradients are calculated. This provides for 


Peete ficam: reduction in the number of £Umction evaluations. 
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